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This  paper  studies  the  thermodynamic  aspects  of  direct  ammonia-fed  solid  oxide  fuel  cells  (DA-SOFC). 
In  this  regard,  anode-supported  oxygen  ion-conducting  solid  oxide  fuel  cells  (SOFC-O)  and  proton¬ 
conducting  solid  oxide  fuel  cells  (SOFC-H)  fuelled  with  ammonia  are  compared  on  the  basis  of  their 
thermodynamic  and  electrochemical  performances  at  open  and  closed  circuit  conditions.  Based  on  the 
thermodynamic  equilibrium,  full  conversion  of  ammonia  is  attained  at  a  temperature  of  700  K  and  atmo¬ 
spheric  pressure.  From  the  results  of  this  study,  it  is  found  that  the  overall  performance  of  the  fuel  cell 
is  largely  dependent  on  critical  factors  such  as  the  operating  temperature  and  the  properties  of  the  elec¬ 
trolyte  and  electrodes  layers.  On  average,  the  peak  power  density  of  ammonia-fed  SOFC-H  is  20-30% 
higher  than  its  SOFC-O  counterpart  under  most  examined  conditions.  This  is  mainly  attributed  to  the 
higher  concentration  of  hydrogen  at  the  anode  of  SOFC-H  in  all  cases. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Trade  market  fluctuations  and  periodic  economic  cycles  like 
the  recent  global  economic  recession  of  2008  have  a  significant 
near-term  impact  on  the  world  energy  demand.  From  a  long-term 
perspective,  there  is  no  doubt  that  the  continuous  growth  of  the 
world  population  presents  great  challenges  to  energy  resources  and 
development.  Therefore,  it  is  evident  that  there  is  an  ever  growing 
need  for  sustainable  and  environmentally  benign  energy  supply  as 
well  as  efficient  power  production  and  distribution. 

Fuel  cells  have  been  identified  as  a  promising  technology  of 
power  production  for  stationary  and  mobile  applications  due 
to  their  high  efficiency  and  small  environmental  footprint.  As 
such,  substantial  research  and  development  efforts  are  invested 
in  advancing  this  technology.  In  particular,  solid  oxide  fuel  cells 
remain  at  the  forefront  of  this  focus  due  to  several  advantages 
like  the  capability  of  directly  operating  on  renewable  and  non¬ 
renewable  energy  carriers  and  fuels. 

As  the  reserve  of  non-renewable  resources  depletes  over  time, 
hydrogen  is  expected  to  emerge  and  attain  a  major  market  share 
as  a  renewable  energy  carrier.  The  shift  of  energy  paradigm  initi¬ 
ated  by  the  new  hydrogen  economy  is  expected  to  have  a  dramatic 
impact  on  the  environmental  and  economic  sustainability  of  soci¬ 
eties  around  the  world.  However,  this  progress  is  hindered  by 
technical  and  economic  challenges  that  must  be  overcome  before 
the  transition  is  realized.  To  that  extent,  research  communities  have 
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started  investigating  the  use  of  suitable  intermediate  hydrogen  car¬ 
riers  in  effort  to  eliminate  these  challenges  and  shortcomings. 

The  task  of  identifying  the  optimal  hydrogen  carrier  is  not 
trivial  as  it  involves  multi-criteria  decision  making  and  consid¬ 
eration  of  various  technical  and  economic  aspects  like  safety, 
energy  density  and  cost  of  processing  or  recycling.  This  has  led 
to  the  careful  scrutiny  of  diverse  spectrum  of  storage  materials 
like  metal  hydrides,  metal-organic  framework  and  amide  systems. 
Despite  extensive  research  and  development  efforts,  these  mate¬ 
rials  and  compounds  have  major  drawbacks  revolving  around  the 
rate  of  hydrogen  desorption,  recyclability  and  high  cost  [  1 ,2].  In  this 
regard,  anhydrous  ammonia  can  potentially  eliminate  the  need  for 
complex  and  perilous  hydrogen  storage  systems. 

Limited  attention  has  been  given  to  the  analysis  of  ammonia- 
fed  solid  oxide  fuel  cells.  For  instance,  Ni  et  al.  [3]  compared 
the  thermodynamic  performance  of  ammonia-fed  oxygen  ion¬ 
conducting  (SOFC-O)  and  proton-conducting  (SOFC-H)  solid  oxide 
fuel  cells  operating  under  temperatures  between  873 1<  and  1273  K. 
The  study  also  assumed  an  indirect  reaction  mechanism  whereby 
ammonia  was  first  decomposed  to  nitrogen  and  hydrogen  while 
the  latter  was  oxidized  to  produce  electric  power.  Comparison  of 
the  cell  potential  of  both  types  of  fuel  cells  over  a  wide  range  of  fuel 
and  oxidant  utilizations  signified  the  superiority  of  SOFC-H  mainly 
due  to  the  absence  of  water  vapour  formation  and  higher  partial 
pressure  of  hydrogen  at  the  anode.  The  maximum  theoretical  (elec¬ 
trochemical)  efficiency  was  found  to  increase  with  increasing  fuel 
utilization  which  is  consistent  with  thermodynamic  studies  on  the 
hydrogen-fed  solid  oxide  fuel  cells  as  given,  for  example,  by  Demin 
and  Tsiakaras  [4]. 

Ni  et  al.  [5]  also  used  an  electrochemical  macro-model  to  com¬ 
pare  the  operating  voltage  and  power  density  of  ammonia-fed 
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Nomenclature 

Bo 

permeability  [m2] 

dx 

depth  of  electrode  layer  [m] 

Deff 

effective  binary  diffusion  [m2  s-1  ] 

Df 

effective  Knudsen  diffusion  [m2  s-1  ] 

E 

reversible  cell  potential  or  voltage  [V] 

E° 

standard  cell  potential  or  voltage  [V] 

/ 

fugacity  of  chemical  species  in  the  system 

f 

fugacity  of  chemical  species  at  STP 

F 

Faraday’s  constant  (96485  C  mol-1 ) 

g 

specific  molar  Gibbs  energy  [J  mol-1  ] 

G 

Gibbs  energy  [J] 

t 

specific  molar  Gibbs  energy  at  STP  [J  mol-1  ] 

j 

current  density  [Am-2] 

Jo 

exchange  current  density  [Am-2] 

M 

molar  mass  [g  mol-1  ] 

n 

number  of  moles  [mol] 

N 

molar  flux  [mol  m-2  s-1  ] 

P 

partial  pressure  [kPa] 

P 

pressure,  total  pressure  [kPa] 

rP 

mean  pore  radius  of  electrode  [m] 

R 

gas  constant  [J  mol-1 1<-1  ] 

Rq 

electrolyte  resistivity  [Q  m] 

T 

temperature  [K] 

U 

utilization  [%] 

V 

working  potential  or  voltage  [V] 

y 

molar  fraction  of  chemical  species 

z 

number  of  electrons  transferred  per  mole  of  fuel 

Greek  letters 

a 

charge  transfer  coefficient 

r 

dimensionless  temperature 

8 

electrolyte  thickness  [m] 

A 

net  change  of  quantity 

£ 

porosity  of  electrode 

e\k 

Lennard -Jones  temperature  parameter  [K] 

Ji 

viscosity  [N  s  m-2] 

G 

molecular  collision  diameter  [A] 

X 

tortuosity  of  electrode 

(Pact, an 

anode  activation  over-potential  [V] 

(Pact,ca 

cathode  activation  over-potential  [V] 

(Pconc.an 

anode  concentration  over-potential  [V] 

( Pconc,ca 

cathode  concentration  over-potential  [V] 

(Pn 

ohmic  over-potential  of  electrolyte  [V] 

collision  integral  of  diffusion 

Subscripts 

eff 

effective 

f 

fuel 

U 

chemical  species 

k 

Knudsen  diffusion 

0 

oxidant 

T 

temperature 

P 

pressure 

Acronyms 

BCNO 

Barium  cerate  neodymium  oxide 

BCSm 

Barium  cerate  samarium  oxide 

DA-SOFC  Direct  ammonia  solid  oxide  fuel  cell 

DeNOx 

Denitrification  (nitrous  oxides) 

DGM 

Dusty  gas  model 

GEM 

Gibbs  energy  minimization 

LHV 

Lower  heating  value 

LSM-YSZ 

;  Lanthanum  strontium  manganite/yttria-stabilized 
zirconia 

Ni-YSZ 

Nickel/yttria-stabilized  zirconia 

OCV 

Open  cell  voltage 

SOFC-H 

Hydrogen  proton  (H+)  conducting  solid  oxide  fuel 
cell 

SOFC-O 

Oxygen  ion  (O2-)  conducting  solid  oxide  fuel  cell 

STP 

Standard  temperature  and  pressure 

TPB 

Triple-phase  boundary 

YSZ 

Yttria-stabilized  zirconia 

SOFC-O  and  SOFC-H  at  a  temperature  of  1073  K  and  atmospheric 
pressure  only.  The  anode-supported  fuel  cells  were  based  on 
yttria-stabilized  zirconia  (YSZ)  and  barium  cerate  samarium  oxide 
(BCSm)  electrolytes  for  SOFC-O  and  SOFC-H,  respectively.  The 
analysis  assumed  a  leak-free  construction,  negligible  interfacial 
resistance,  no  fuel  cross-over  or  electronic  conductivity  through 
the  electrolyte  layer.  Thus  only  activation,  ohmic  and  concentration 
over-potentials  were  addressed.  Their  work  also  emphasized  the 
higher  performance  of  SOFC-O  over  SOFC-H.  This  was  explained 
by  the  high  ohmic  over-potential  of  BCSm  electrolyte  and  the  high 
concentration  over-potential  at  the  SOFC-H  cathode  as  the  water 
vapour  formation  obstructed  the  oxygen  transport  to  the  reaction 
sites.  The  results  were  quite  different  from  those  presented  by 
Assabumrungrat  et  al.  [6]  and  Biesheuvel  and  Geerlings  [7]  who 
advocated  the  superior  performance  of  SOFC-H  in  comparison  to 
SOFC-O  using  methanol  and  methane  fuels  at  the  same  operating 
conditions. 

With  the  focus  on  SOFC-H,  another  study  was  conducted  by  Ni 
et  al.  [8]  to  examine  the  impact  of  different  types  of  structural  sup¬ 
port  on  the  performance  of  the  fuel  cell.  It  was  found  that  higher 
power  density  can  be  achieved  by  the  anode-supported  cell  pri¬ 
marily  because  of  the  thinner  electrolyte  material  which  greatly 
reduced  the  ohmic  polarization.  Later,  their  electrochemical  model 
was  improved  [9]  by  incorporating  a  first-order  model  of  the  kinetic 
rate  of  ammonia  decomposition  which  was  previously  developed 
by  Chellappa  et  al.  [10].  The  kinetic  model  was  coupled  with  the 
dusty  gas  model  (DGM)  to  obtain  the  molar  fraction  of  ammo¬ 
nia  across  the  anode  thickness  at  various  temperatures.  The  model 
showed  that  the  ammonia  was  fully  decomposed  at  an  approximate 
anode  depth  of  25  p,m  and  300  p>m  at  1073  K  and  873  K,  respec¬ 
tively.  It  was  also  shown  that  ammonia  did  not  fully  decompose  at 
temperatures  below  873 1<  which  resulted  in  a  low-partial  pressure 
of  hydrogen  and  deteriorated  performance  of  the  fuel  cell. 

This  study  is  motivated  by  the  technical  challenges  to  the 
rapid  development  of  hydrogen  markets  and  infrastructure  which 
include  the  constraints  of  safe  transportation,  low  temperature 
storage  as  well  as  the  high  pressure  necessary  to  achieve  the 
required  energy  density.  The  objective  of  this  paper  is  to  high¬ 
light  the  advantages  and  merits  of  ammonia  as  an  energy  carrier 
by  shedding  light  on  various  important  aspects  like  its  energy  stor¬ 
age  capacity  as  well  as  its  demand  and  supply  around  the  world. 
Furthermore,  the  aim  of  this  work  is  to  examine  the  operating 
principles  of  ammonia-fed  solid  oxide  fuel  cells  and  investigate 
the  cell-level  thermodynamic  and  electrochemical  performances 
at  open  and  closed  circuit  conditions. 

2.  Ammonia  as  a  hydrogen  carrier 

Ammonia  is  extensively  used  in  many  industries  including 
refrigeration,  petrochemical  and  food  processing.  Furthermore, 
large  quantities  of  ammonia  are  used  worldwide  for  agricultural 
purposes.  It  is  estimated  that  the  global  capacity  of  ammonia 


Table  1 

Comparison  of  energy  storage  capacity. 
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Energy  carrier 

Density  (gL-1) 

Fraction  of  hydrogen  (%  H2) 

Volumetric  density  (gH2  L-1 ) 

Energy  density  (MJ  L-1 ) 

Gaseous  H2  (298  K,  10  MPa) 

7.68 

100 

7.7 

0.92 

Liquid  H2  (20  K,  0.1  MPa) 

71.28 

100 

71.3 

8.56 

Liquid  NH3  (298  K,  1  MPa) 

603 

17.76 

107.1 

12.85 

Sources:  [15,16]. 

Energy  density  is  based  on  LHV  of  hydrogen. 


production  was  about  153  million  tonnes  in  2009  and  is  expected 
to  increase  to  244  million  tonnes  by  2014  [11].  This  major  increase 
is,  in  part,  due  to  the  growing  use  of  ammonia  in  emerging  emission 
control  (DeNOx)  technologies  in  industrial  and  automotive  appli¬ 
cations. 

The  infrastructure  and  technology  of  ammonia  production  are 
also  well  established  with  existing  industrial  plants  around  the 
world  to  support  the  increasing  demand  for  ammonia.  Steam 
reforming  of  natural  gas  is  the  current  primary  method  of  ammo¬ 
nia  production.  Consequently,  the  best  available  technologies  emits 
about  1.7  tonnes  of  carbon  dioxide  per  tonne  of  ammonia  produc¬ 
tion  most  of  which  can  be  easily  recovered  for  use  in  downstream 
processes  such  as  the  production  of  urea  or  other  derivatives  [12]. 
This  figure  excludes  the  potential  amount  of  carbon  dioxide  emit¬ 
ted  when  carbon-based  fuel  is  used  to  provide  the  motive  power 
to  the  process  of  ammonia  production. 

Recently,  ammonia  has  been  heralded  as  an  excellent  hydrogen 
carrier  for  its  several  favourable  attributes.  For  example,  it  can  be 
liquefied  at  ambient  temperature  under  a  moderate  pressure  of 
1.0  MPa  for  the  purpose  of  storage  and  transportation.  As  shown 
in  Table  1,  its  energy  density  is  48%  higher  than  that  of  liquefied 
hydrogen.  The  gravimetric  and  volumetric  capacities  of  ammonia 
also  exceed  the  ultimate  target  for  hydrogen  storage  capacity  set  by 
the  US  Department  of  Energy  [13].  Furthermore,  ammonia  leakage 
is  easily  detectable  at  concentrations  as  low  as  17  ppm  [14]  due  to 
its  pungent  odour. 

In  contrast,  the  utilization  of  ammonia  in  certain  applications 
is  hindered  by  safety  implications  among  other  concerns.  Ammo¬ 
nia  vapour  can  be  toxic  and  life-threatening  at  high  concentration. 
It  is  also  known  that  exposure  to  the  vapour  at  a  concentration  of 
1700  ppm  can  cause  severe  lung  damage  or  respiratory  failure  and 
death  at  5000  ppm  [  1 7 ].  Another  significant  concern  is  the  corrosive 
effects  of  ammonia  on  some  materials  like  copper,  zinc  and  their 
alloys.  In  order  to  mitigate  these  risks  and  shortcomings,  some 
focus  has  been  directed  towards  stabilizing  the  ammonia  by 
combining  it  in  metal  ammine  complexes  [18].  This  allows  for 
transportation  and  long-term  storage  of  ammonia  in  solid-state 
form  and  provides  the  means  of  releasing  the  hydrogen  as  needed. 


3.  Analysis 

3.1.  Thermodynamic  equilibrium  of  ammonia  decomposition 

When  heated,  ammonia  decomposes  to  hydrogen  and  nitrogen 
according  to  the  following  reaction: 

nh3^|h2  +  1n2  (i) 

The  above  reaction  is  endothermic  and  depends  on  the  temper¬ 
ature  and  pressure  of  the  system  among  other  factors.  The  reaction 
proceeds  to  the  right  as  the  temperature  increases  or  to  the  left  as 
the  pressure  increases  thus  producing  less  hydrogen.  The  thermo¬ 
dynamic  equilibrium  of  ammonia  can  be  obtained  by  minimizing 
the  total  Gibbs  energy  of  the  reaction: 


(AGsystem  )jp  —  0 


(2) 


where  the  Gibbs  energy  of  the  system  can  be  defined  as  the  sum¬ 
mation  of  the  product  of  moles  of  chemical  species  i  multiplied  by 
the  respective  specific  Gibbs  energy  at  constant  temperature  and 
pressure 

^system  —  Zni&i  (3) 

The  specific  Gibbs  energy  can  be  described  as  the  sum  of  the 
standard  Gibbs  function  of  formation  and  the  chemical  activity  of  a 
given  species  i.  For  real  gases,  the  activity  of  species  is  defined  as  the 
ratio  of  its  fugacity  in  the  system  to  that  at  standard  temperature 
and  pressure  (STP)  conditions: 

gi  =gfi  +  RT  ln^  (4) 

•U 

If  ideal  or  perfect  gas  conditions  are  assumed,  the  activity  of  the 
gaseous  species  is  equal  to  its  partial  pressure  in  the  system.  For 
condensed  matter  (solid  or  liquid),  the  activity  is  equal  to  unity. 
Therefore,  the  total  Gibbs  energy  of  the  system  is 

Csystem  =  (V  n,[g£  +  RT  ln(y,P)])  +  (5) 

'  /  gas  \z — '  /  condensed 

As  an  advantage  of  this  method,  chemical  reactions  need  not  to 
be  specified.  Instead,  Gibbs  energy  minimization  (GEM)  requires 
data  of  the  standard  Gibbs  function  of  formation  of  all  chemical 
species  which  are  expected  to  be  involved  or  contribute  to  the 
thermodynamic  equilibrium  under  consideration. 

3.2.  Description  and  operational  principles  ofDA-SOFC 

The  operating  principle  of  direct  ammonia-fed  solid  oxide  fuel 
cell  (DA-SOFC)  with  oxygen  ion  conducting  electrolyte  is  assumed 
to  take  an  indirect  path  to  the  oxidation  of  the  fuel  by  thermally 
decomposing  the  ammonia  to  extract  the  hydrogen  which  will  be 
electrochemically  oxidized  in  the  fuel  cell.  Depending  on  the  ther¬ 
modynamic  equilibrium  of  the  ammonia  decomposition  reaction  at 
given  temperature  and  pressure,  nitrogen  gas  may  be  produced  in 
high  concentration.  The  inert  nitrogen  gas  can  dilute  the  hydrogen 
concentration  at  the  anode  resulting  in  reduction  of  the  reversible 
cell  potential.  Furthermore,  the  nitrogen  gas  produced  is  discharged 
from  the  anode  compartment  along  with  other  gaseous  effluents  as 
depicted  in  Fig.  1. 

On  the  cathode  side,  oxygen  from  air  is  transported  across 
the  cathode  layer  to  the  cathode-electrolyte  interface  where  it 
is  reduced  to  oxygen  ion.  The  oxygen  ions  are  then  transported 
through  the  electrolyte  to  electrochemically  react  with  hydrogen 


Fig.  1.  Schematic  representation  of  oxygen  ion-conducting  direct  ammonia  solid 
oxide  fuel  cell. 
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NH3  ■= 

Anode 

Electrolyte 

Cathode 

02,  N2  i= 


=^>  NH3  ->  3/2  H2  +  1/2  N2  c 

==>  H 2, 

h2*  |h2 

Ih* 

IrT- 

o2f  |h2o  h2o| 

to2 

=>  02 

Fig.  2.  Schematic  representation  of  proton-conducting  direct  ammonia  solid  oxide 
fuel  cell. 


at  the  anode-electrolyte  interface.  The  water  vapour  produced  by 
the  reaction  is  transported  to  the  anode  channel  which  can  also 
dilute  the  fuel  concentration.  The  water  vapour  then  exits  the  cell 
along  with  unreacted  fuel  and  other  effluents. 

Fig.  2  depicts  the  configuration  of  direct  ammonia-fed  solid 
oxide  fuel  cell  (DA-SOFC)  with  proton-conducting  electrolyte.  As 
ammonia  enters  the  cell  at  the  anode  side,  it  decomposes  to 
hydrogen  and  nitrogen  in  proportions  governed  by  the  thermo¬ 
dynamic  equilibrium  at  specific  temperature  and  pressure.  While 
nitrogen  acts  as  an  inert  diluent,  hydrogen  is  transported  to 
the  anode-electrolyte  interface  to  be  oxidized  to  hydrogen  pro¬ 
tons.  Subsequently,  hydrogen  protons  are  transported  through  the 
electrolyte  layer  to  the  cathode-electrolyte  interface  to  undergo 
electrochemical  reaction  with  oxygen  from  the  cathode  side.  The 
water  vapour  produced  by  the  fuel  cell  reaction  is  then  transferred 
to  the  bulk  gas  in  the  cathode  channel  where  it  exits  the  cell  along 
with  other  effluents.  This  is  a  major  advantage  since  hydrogen  is  not 
diluted  by  the  water  vapour  generated  from  the  fuel  cell  reaction. 


3.3.  Thermodynamic  and  electrochemical  analyses 


The  equilibrium  potential  of  any  fuel  cell  is  affected  by  the  partial 
pressures  of  the  species  involved  in  the  electrochemical  reaction. 
Under  open-circuit  conditions,  the  potential  (voltage)  can  be  writ¬ 
ten  as 


E  =  E 


,  RT. 

+  2 Fln 


Ph2(Po2)1/2 

Ph2o 


(6) 


where 
o  -AG 

E  =^r 


(7) 


Furthermore,  the  examination  of  the  potential  of  solid  oxide  fuel 
cells  at  closed  circuit  conditions  requires  the  characterization  of  the 
over-potentials  affecting  the  cell  operation.  The  working  potential 
or  voltage  of  SOFC  can  be  characterized  as 


V  —  E  —  (Pact,  an  ~  (Pact,ca  ~  ~  cone,  an  ~  (Pconc,ca 


(8) 


The  activation  over-potential  is  an  irreversible  loss  related  to  the 
electrode  kinetics  and  represents  the  difference  between  the  actual 
potential  (closed-circuit  conditions)  and  the  reversible  potential 
(open-circuit  conditions).  The  implicit  relationship  between  the 
current  density  and  the  activation  over-potential  can  be  described 
using  the  Butler-Volmer  expression  [19] 


J=Jo 


exp 


/ oczF(pact\ 
V  RT  ) 
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f  (1  -  a)zF(pact 
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(9a) 


The  exchange  current  density  (J0)  is  a  measure  of  the  magnitude 
of  electron  activity  at  the  equilibrium  potential  of  the  electrode.  Its 
value  is  highly  influenced  by  the  electrode  material,  structure  and 
other  factors  such  as  the  temperature  of  the  reaction  and  the  length 
of  the  triple-phase  boundary  (TPB)  [20].  The  charge  transfer  coef¬ 
ficient  (a)  describes  the  effect  of  the  electric  potential  on  the  ratio 
of  the  forward  to  the  reverse  activation  barrier.  For  most  electro¬ 
chemical  reactions  in  solid  oxide  fuel  cells,  this  value  is  assumed 


to  be  0.5  [21].  As  such,  the  explicit  relationship  of  the  activation 
over-potential  is 

R.T  (  I  \ 

(pact  y  —  -pr  sin  ft-1  (  -y—  )  ,  y  =  anode,  cathode  (9b) 

b  \zJ0,y  J 

The  electronic  resistance  of  the  fuel  cell  electrode  and  inter¬ 
connect  is  usually  very  small  and  is  considered  negligible  when 
compared  to  the  ionic  resistance  of  the  electrolyte.  The  ionic  resis¬ 
tance  can  be  described  by  Ohm’s  law  [22]  as 

<Pa=Mn  (10) 

The  electrolyte  resistance  (ft^)  is  a  function  of  the  electrolyte 
properties  and  is  highly  dependent  on  the  operating  temperature. 

The  concentration  over-potential  accounts  for  losses  incurred  by 
the  resistance  of  the  porous  electrode  to  the  transport  of  gaseous 
species  between  the  gas  channel  and  the  reaction  sites  at  TPB.  The 
mass  transport  in  the  electrodes  is  driven  by  the  diffusion  of  react¬ 
ing  species  due  to  concentration  gradient  as  well  as  the  permeation 
caused  by  the  pressure  gradient.  Several  mass  transport  models 
with  varied  accuracies  have  been  used  to  capture  the  effects  of 
concentration  over-potential  on  the  performance  of  fuel  cells. 

It  has  been  demonstrated  that  the  highest  accuracy  can  be 
achieved  using  the  dusty  gas  model  [23]  and  the  mean  transport 
pore  model  (MTPM)  [24].  In  particular,  DGM  visualizes  the  porous 
medium  as  a  collection  of  large  spherical  particles  suspended  in 
space  by  external  forces.  The  model  takes  into  account  the  effects 
of  diffusion  and  permeation  mechanisms  on  the  concentration  of 
reacting  species  which  is  essential  when  modelling  pressurized  fuel 
cell  systems.  The  one-dimensional  multi-component  mass  trans¬ 
port  DGM  can  be  written  as 
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where  the  effective  Knudsen  diffusion  accounts  for  the  porosity  and 
tortuosity  of  the  electrode  such  that  [25]: 


eff  =  2  £  /8RT 

L  k  3  t  y  jrMj  p 


(12) 


The  effective  binary  diffusion  of  chemical  species  can  be 
obtained  using  the  Chapman-Enskog  equation  [26]: 


Df 


0.0018583- 

r 


(h-  1  Mj)  P<T?.C2DJj 


(13) 


and  the  mean  characteristic  length  (cry)  of  the  molecular  collision 
diameter  of  species  i  and  j  is  given  by 


fa  +  cry) 
2 


(14) 


The  collision  integral  (T2Dij)  is  a  function  of  temperature  and  the 
Lennard -Jones  parameter  (k/Sy)  such  that 

1.06036  0.19300  1.03587 

DJj  ~  r°  i5610  +  exp(0.47635E)  +  exp(  1.52996/") 
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The  values  of  the  parameters  in  Eqs.  (12)-(17)  have  been 
obtained  from  Refs.  [26,27].  The  permeability  of  the  porous  elec¬ 
trode  is  estimated  using  the  Kozeny-Carman  relationship  [28] 


4s3  r3 

72t(1  -  e)2 


(16') 


Furthermore,  the  viscosity  of  gas  mixture  can  be  calculated  using 
the  following  semi-empirical  equation  [29]: 


P'  mix  — 


E  yty/Mi 


(17') 


The  gas  viscosities  are  obtained  from  Ref.  [26].  As  presented  by 
Zhu  et  al.  [30],  the  summation  of  Eq.  (1 1 )  allows  for  the  evaluation 
of  the  pressure  gradient  (dP/dx)  across  the  electrode  layer  as 


dp  _  Ej^W/gg) 

dx  (1  /RT)  +  (B0P/nRT)J2i=1(yi/Dfk) 


At  TPB,  the  molar  flux  of  the  gaseous  reactants  involved  in  the 
electrochemical  reaction,  namely  hydrogen  and  oxygen,  can  be 
related  to  the  current  density  as 

Ni  =  If  (19) 


Depending  on  the  type  of  electrolyte  (SOFC-O  or  SOFC-H),  the 
molar  flux  of  water  vapour  can  be  calculated  using  Graham’s  law  of 
diffusion  which  governs  the  diffusion  process  of  gas  mixtures  [23]: 


Y^Niy/Ml  =  0 


(20) 


The  molar  fluxes  of  all  other  non-reacting  species  are  equal  to 
zero  [5].  Finally,  the  relationship  between  the  concentration  over¬ 
potential  and  the  partial  pressures  can  be  written  as  [31  ] 
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Fig.  3.  Thermodynamic  equilibrium  of  ammonia. 


4.  Results  and  discussion 

In  order  to  obtain  the  molar  fractions  and  partial  pressures  of  the 
reacting  species  at  TPB,  a  computer  code  has  been  developed  using 
function  (ode45)  in  MATLAB  (MathWorks,  Natick,  MA,  USA)  [32] 
to  numerically  solve  the  set  of  simultaneous  differential  equations 
represented  by  Eq.  (11).  The  required  thermodynamic  data  were 
obtained  from  the  JANAF  tables  [33].  The  minimization  of  Eq.  (5) 
was  performed  using  the  optimization  function  in  MATLAB. 


4.1.  Thermolysis  of  ammonia  in  DA-SOFC 

The  use  of  the  GEM  method  revealed  that  the  decomposition 
of  ammonia  is  favoured  at  high  temperature  and  low  pressure 
as  shown  in  Fig.  3.  In  particular,  full  conversion  of  ammonia  can 
be  achieved  at  a  temperature  of  700  K  and  atmospheric  pressure. 
Nonetheless,  it  is  well  known  that  the  kinetic  rate  of  reaction  is 
rather  slow  and  often  requires  the  use  of  a  suitable  catalyst  to 
promote  efficient  conversion  of  ammonia.  Complete  conversion 
may  be  achieved  over  expensive  metal  catalysts  such  as  ruthenium 
and  platinum  at  temperatures  around  600-650  K  [34];  however, 
more  economic  alternatives  include  iron  and  nickel-based  cata¬ 
lysts.  Since  the  operating  temperature  of  solid  oxide  fuel  cells 
commonly  ranges  between  900  and  1400  K  [35],  it  is  safe  to  neglect 
the  kinetic  rate  of  reaction  and  assume  that  full  decomposition  of 
ammonia  is  attained  within  the  porous  anode  layer  of  the  fuel  cell. 


4.2.  Reversible  cell  potential  (open  cell  voltage) 


3.4.  Fuel  and  oxidant  utilization 


The  utilization  of  ammonia  or  urea  can  be  defined  in  terms 
of  the  actual  supply  and  consumption  of  the  fuel  or  its  hydrogen 
equivalent  such  that 


Uf  = 


(Fuel)consumed 

(Fuel)supplied 


(H2  Consumed 

(^2  Supplied 


(22) 


Similarly,  the  utilization  of  air  as  an  oxidant  can  be  expressed 
as  a  function  of  the  air  supply  and  usage  or  its  oxygen  equivalent 
since  air  can  be  assumed  to  be  composed  of  79%  nitrogen  and  21% 
oxygen  by  volume. 


U0  = 
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(23) 


Fig.  4  illustrates  the  variation  of  the  molar  fraction  of  hydro¬ 
gen  as  a  function  of  ammonia  utilization  at  temperature  of  1073  K 
and  atmospheric  pressure.  At  the  anode  of  SOFC-O,  hydrogen 
is  balanced  by  the  equimolar  formation  of  water  vapour  due  to 
the  electrochemical  reaction  which  results  in  a  linear  profile  of 
the  hydrogen  molar  fraction.  For  SOFC-O,  a  non-linear  profile  is 
observed  mainly  due  to  the  decomposition  of  ammonia  and  con¬ 
sumption  of  hydrogen. 

Fig.  5  shows  that  the  molar  fraction  of  water  vapour  at  the  cath¬ 
ode  of  SOFC-O  is  higher  than  that  at  the  anode  of  SOFC-H  at  any 
ammonia  utilization  factor  above  13%.  This  profile  is  observed  due 
to  the  production  of  water  vapour  in  quantity  equal  to  the  hydro¬ 
gen  consumption  at  the  anode  of  SOFC-O.  Conversely,  the  constant 
molar  fraction  of  water  vapour  formation  at  the  cathode  of  SOFC-H 
is  dependent  solely  on  the  invariant  oxidant  utilization. 
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It  is  always  desirable  to  have  a  high  fuel  utilization  to  maximize 
the  system  efficiency;  however,  high  open  cell  voltage  (OCV)  is  a 
competing  objective  which  tends  to  decline  as  the  fuel  utilization  is 
increased.  Fig.  6  exhibits  the  OCV  of  SOFC-O  and  SOFC-FI  as  a  func¬ 
tion  of  fuel  utilization.  It  can  be  seen  that  the  difference  becomes 
more  pronounced  at  higher  fuel  utilizations  due  to  the  depletion  of 
hydrogen  and  formation  of  water  vapour. 

It  is  also  worth  noting  that  the  prediction  of  OCV  values  in  this 
study  is  consistent  with  the  experimental  findings  of  Fuerte  et  al. 
[36]  who  measured  the  voltage  of  SOFC-O  fuelled  with  hydrogen 
and  ammonia  at  temperatures  of  973-1 1 73  K.  Fig.  7  shows  the  OCV 
degradation  of  both  hydrogen  and  ammonia-fed  fuel  cells  as  the 
operating  temperature  of  the  cell  increases  which  also  indicates 
that  the  reaction  mechanisms  of  both  cells  are  similar  and  limited 
to  the  direct  electrochemical  oxidation  of  hydrogen. 

Fig.  8  indicates  that  the  OCV  can  be  improved  by  5-6%  when 
the  operating  pressure  is  increased  by  five  times  the  atmospheric 
pressure.  This  modest  enhancement  is  of  particular  benefit  when 
the  fuel  cell  is  integrated  with  gas  turbine  systems  (SOFC/GT). 

In  addition,  increasing  the  air  utilization  reduces  the  OCV  as 
shown  in  Fig.  9.  The  deterioration  is  more  pronounced  in  SOFC-FI 
due  to  higher  partial  pressure  of  water  vapour  and  lower  partial 
pressure  of  oxygen  at  the  cathode. 


1.0  n 

T=  1073  K  - SOFC-H 

P=  101  kPa  - SOFC-O 

Air  Utilization  =  25% 

0.8- 


Ammonia  Utilization  [%] 


l.l 


1.0- 


0.7- 


SOFC-H 
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Ammonia  Utilization  [%] 


Fig.  6.  Open  cell  voltage  of  DA-SOFC  at  different  temperatures. 
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Fig.  5.  Molar  fraction  of  water  vapour  in  DA-SOFC. 


Fig.  8.  Effect  of  pressure  on  the  open  cell  voltage  of  DA-SOFC. 
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Fig.  9.  Effect  of  oxidant  (air)  utilization  on  the  open  cell  voltage  of  DA-SOFC. 


Fig.  10.  Comparison  of  simulation  and  experimental  results  of  ammonia-fed 
SOFC-O. 


4.3.  Operational  voltage  and  performance 
4.3 A.  Model  validation 

In  the  open  literature,  experimental  investigations  on  ammonia- 
fed  solid  oxide  fuel  cells  are  scarce  and  rarely  contain  all  the 
necessary  information  and  parameters  needed  for  model  evalua¬ 
tion.  Therefore,  certain  assumptions  were  made  and  some  model 
parameters  were  adjusted  to  agree  with  the  experimental  data  as 
reflected  in  Table  2.  Ma  et  al.  [  1 5  ]  conducted  an  experiment  in  which 
ammonia  was  used  as  a  fuel  in  anode-supported  SOFC-O  with  the 
arrangement  of  Ni-YSZ|YSZ|LSM-YSZ.  The  cell  was  tested  at  atmo¬ 
spheric  pressure  and  different  temperatures  between  923  K  and 
1123  K.  It  is  worth  noting  that  experimental  tests  are  typically  per¬ 
formed  with  low  utilization  of  reactants  in  order  to  avoid  the  risk  of 
local  overheating  and  degradation  which  can  occur  at  high  current 
densities.  In  practice,  higher  utilization  factors  are  always  desirable 
despite  the  lower  efficiency  and  specific  power  output  of  the  fuel 
cell  stack. 

In  addition,  the  performance  of  the  cell  was  characterized  by 
the  high  interfacial  resistance  which  accounted  for  about  50%  of 
the  total  resistance  of  the  cell.  The  resistance  of  the  electrolyte  was 
also  higher  than  expected  at  the  same  temperature  as  reported 
elsewhere  [31,37].  In  part,  this  can  be  attributed  to  the  poor 


Table  2 

Parameters  used  in  the  validation  of  ion-conducting  DA-SOFC. 


Parameters 

Values 

Operating  temperature,  T  (K) 

1123 

Operating  pressure,  P  (kPa) 

101 

Fuel  utilization,  Uf  (%) 

5 

Oxidant  (air)  utilization,  U0  (%) 

5 

Exchange  current  density  of 
anode  Jan  (Am-2) 

7x1°9(>)  (^)exp(-^) 

Exchange  current  density  of 

7xl0*(>y'25exppr°) 

cathode Jca  (Am-2) 

Electrode  porosity,  s 

0.4 

Electrode  tortuosity,  r 

4.25 

Average  pore  radius  of 

0.5 

electrode,  rp  (p,m) 

Resistivity  of  electrolyte, 

0.735 

(£2m) 

Anode  thickness  (|jim) 

500 

Electrolyte  thickness  (p,m) 

30 

Cathode  thickness  (|jim) 

10 

contact  between  the  electrodes  and  the  electrolyte  as  well  as  the 
microstructural  properties  of  the  cell  components  which  are  highly 
dependent  on  the  processing  method  as  discussed  by  Ringuede 
et  al.  [38].  The  results  in  Fig.  10  show  excellent  agreement  (within 
2%)  between  the  numerical  simulation  and  the  experimental  data. 

4.3.2.  Comparison  of  ion  and  proton-conducting  DA-SOFC 

Certain  modelling  parameters  have  been  fixed  to  establish  a 
common  basis  of  comparison.  The  parameters  for  oxygen  ion  and 
proton-conducting  DA-SOFC  are  listed  in  Table  3.  Fig.  1 1  shows  the 
polarization  curve  of  DA-SOFC  operating  at  1 073  K  and  atmospheric 
pressure.  It  can  be  seen  that  high  current  density  can  be  achieved 
by  employing  an  anode-supported  arrangement  of  the  fuel  cell. 
This  can  effectively  decrease  the  ohmic  over-potential  due  to  the 
reduction  in  the  electrolyte  thickness.  The  peak  power  density  of 
SOFC-H  is  about  21%  higher  than  that  of  SOFC-O.  In  addition,  the 
voltage  decreases  rapidly  for  both  types  of  fuel  cells  at  high  current 
densities  due  to  the  low  concentrations  of  reactants  at  TPB. 


Table  3 

Parameters  used  in  modelling  the  polarization  curve  of  DA-SOFC. 


Parameters 

Values 

Fuel  utilization,  Uf  (%) 

5 

Oxidant  utilization,  U0  (%) 

5 

Exchange  current  density  of  anode 

7  x 

ISOFC— OJan  (Am-2) 

109(^f)  (^)  exp  (-122000) 

Exchange  current  density  of 

7xl°9(^)025exP  (-t) 

cathode  |SOFC-OJca  (Am-2) 

Exchange  current  density  of  anode 

2500@1073I< 

|SOFC-H,Jan  (Am-2) 

5000@1273I< 

Exchange  current  density  of 

1700@1073I< 

cathode  |SOFC-HJca  (Am-2) 

3400@1273I< 

Electrode  porosity,  e 

0.4 

Electrode  tortuosity,  r 

4.25 

Average  pore  radius  of  electrode,  rp 

0.5 

(pan) 

Resistivity  of  electrolyte  |SOFC-0, 

2.94  x  10-5  exp  ^  ^ 

Rq  (£2m) 

Resistivity  of  electrolyte  |SOFC-H, 

0.37@1073I< 

Rq  (£2  m) 

0.21@1273I< 

Anode  thickness  (|jim) 

500 

Electrolyte  thickness  (p,m) 

30 

Cathode  thickness  (|Jim) 

30 

Sources:  [8,15,39,40]. 


Sources:  [8,15,37,39-42]. 
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Current  Density  [A  m'2] 

Fig.  11.  Polarization  curve  of  DA-SOFC  at  1073  K  and  atmospheric  pressure. 

As  depicted  in  Fig.  12,  the  activation  over-potentials  of  SOFC-O 
and  SOFC-H  anodes  are  comparable;  however,  it  is  rather  higher  at 
the  cathode  of  SOFC-FI  in  comparison  to  that  at  the  cathode  of  its 
counterpart.  This  is  highly  affected  by  the  difference  in  exchange 
current  density  which  is  estimated  to  be  about  2200 Am-2  for 
SOFC-O  and  1700  Am-2  for  SOFC-FI  at  the  given  temperature  and 
pressure.  The  effect  of  higher  activation  over-potential  in  SOFC-H 
is  compensated  by  the  low  ohmic  over-potential  as  shown  in 
Fig.  13.  The  resistivity  (or  ionic  conductivity)  of  electrolyte  is  highly 
influenced  by  the  chemical  composition,  processing  method  and 
sintering  temperature  [38].  The  values  of  resistivity  used  in  this 
numerical  study  are  0.455  £2  m  [43  ]  for  oxygen  ion-conducting  elec¬ 
trolyte  (YSZ)  and  0.370  £2m  [39]  for  proton-conducting  barium 
cerate  neodymium  oxide  (BCNO)  electrolyte  at  1073K  and  repre¬ 
sent  an  average  of  the  values  found  in  the  literature. 

Fig.  14  reflects  the  effect  of  concentration  over-potential  on  the 
performance  of  anode-supported  DA-SOFC.  Due  to  the  thickness  of 
the  anode,  its  concentration  over-potential  is  significantly  higher 
than  that  of  the  cathode.  Moreover,  the  over-potential  is  more  pro¬ 
nounced  at  the  anode  of  SOFC-O  due  to  the  formation  of  water 
vapour  at  the  anode-electrolyte  interface  and  its  diffusion  through 
the  porous  anode  into  the  gas  channel. 

Additionally,  the  operating  temperature  has  a  major  impact  on 
the  performance  of  the  fuel  cell.  Fig.  15  demonstrates  the  increase 
in  power  density  despite  the  slight  decrease  in  OCV  caused  by  the 


increase  in  operating  temperature  from  1073  to  1273  K.  The  higher 
operating  voltage  and  current  density  are  the  result  of  lower  over¬ 
potentials  at  high  temperatures.  The  peak  power  density  of  SOFC-O 
is  about  50%  higher  than  that  of  SOFC-H  which  is  mainly  due  to  the 


Fig.  15.  Polarization  curve  of  DA-SOFC  at  1273  K  and  atmospheric  pressure. 
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Fig.  16.  Polarization  curve  of  DA-SOFC  at  1073  K  and  pressure  of  505  kPa. 


uncertainty  in  estimating  the  exchange  current  density  of  SOFC-H 
in  contrast  to  the  rate  expression  used  for  SOFC-O  counterpart. 
The  difference  in  peak  power  density  obtained  in  this  study  may  be 
somewhat  justified  by  recognizing  that  the  performance  of  SOFC-FI 
deteriorates  at  elevated  temperatures  due  to  the  high  electronic 
conductivity  of  the  electrolyte  as  given  by  various  accounts  in  the 
literature  (e.g.  [39,42]).  Therefore,  it  is  highly  discouraged  to  utilize 
and  operate  SOFC-FI  cells  at  temperatures  beyond  1073  K. 

Finally,  Fig.  16  illustrates  the  increase  in  the  power  density 
of  the  fuel  cell  at  elevated  operating  pressure.  The  performance 
improvement  is  the  result  of  a  slight  increase  in  OCV  and  moderate 
enhancement  of  the  permeation  of  gaseous  species  to  and  from  the 
reaction  sites  at  TPB. 

5.  Conclusions 

The  merits  of  using  ammonia  as  a  hydrogen  carrier  have  been 
highlighted  in  this  study.  The  operating  principles  of  ammonia- 
fed  solid  oxide  fuel  cells  have  been  illustrated  as  well.  In  order 
to  investigate  the  cell-level  thermodynamic  and  electrochemical 
performances  at  open  and  closed  circuit  conditions,  a  model  was 
developed  and  validated  by  excellent  agreement  with  experimental 
data  obtained  from  the  literature. 

Furthermore,  the  parametric  analyses  performed  in  this  work 
revealed  that  proton-conducting  DA-SOFC  maintained  higher 
reversible  cell  potential  than  the  ion-conducting  counterpart  under 
all  open-circuit  conditions.  Under  closed-circuit  conditions,  the 
proton-conducting  DA-SOFC  demonstrated  a  better  performance  in 
comparison  to  the  ion-conducting  counterpart  due  to  higher  partial 
pressure  of  hydrogen  at  the  anode  under  all  operating  conditions. 

The  results  of  this  work  emphasized  the  advantages  of  the 
operating  principle  and  performance  of  SOFC-FI;  however,  their 
performance  is  highly  influenced  by  many  factors  such  as  the  type 
of  electrolyte  used  and  the  microstructural  properties  of  the  elec¬ 
trodes.  As  such,  all  design  parameter  must  be  optimized  to  achieve 
the  best  possible  performance. 
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